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We have analyzed the Kuramoto-Sivashinsky equation with a stochastic noise term through a 
dynamic renormalization group calculation. For a system in which the lattice spacing is smaller 
than the typical wavelength of the linear instability occurring in the system, the large-distance and 
long-time behavior of this equation is the same as for the Kardar-Parisi-Zhang equation in one 
and two spatial dimensions. For the d = 2 case the agreement is only qualitative. On the other 
hand, when coarse-graining on larger scales the asymptotic flow depends on the initial values of the 
parameters. 



PACS numbers: 64.60-Ht, 68.35. Rh, 05.40. +j, 79.20.Rf 



I. INTRODUCTION 

Currently, there is much interest in understanding 
the formation and roughening of nonequilibrium inter- 
faces A common feature of many interfaces ob- 
served experimentally and in discrete models is that their 
roughening follows simple scaling laws characterized by 
the roughness exponent a and the dynamic exponent z, 
which determine the scaling behavior of, e.g., the correla- 
tion functions. In many cases, the scaling exponents can 
be obtained using stochastic evolution equations of which 
a seminal example is the so-called Kardar-Parisi-Zhang 
(KPZ) equation 

Another equation, which has been actively discussed 
in problems of pattern formation, such as chemical tur- 
bulence and flame-front propagation, is the so-called 
Kuramoto-Sivashinsky (KS) equation This is 

a deterministic nonlinear equation which exhibits spa- 
tiotemporal chaos. Qualitatively, the chaotic nature of 
the KS equation generates stochasticity in such a way 
that its solution displays scaling at large distances and 
long times. An important question to answer is whether 
the KPZ and KS equations belong to the same or to dif- 
ferent universality classes, i.e., whether the scaling prop- 
erties of the interfaces are described by the same or differ- 
ent values for the critical exponents. In 1 + 1 dimension 
there are numerical ]F^-|ll| and analytical ]l^-|l9|] results 
which show that the KS and KPZ equations indeed ex- 
hibit the same scaling behavior. However, it is still an 
open question whether the KS and the KPZ equations 
fall into the same universality class in 2 + 1 dimensions 

In order to address these questions from a different 
point of view, we have investigated a noisy version of the 
KS equation by a dynamic renormalization group (RG) 
analysis. Even though the relation of such a noisy KS 
equation to the KS system is still to be completely clari- 
fied, we believe that the results we obtain for the scaling 
behavior of the noisy KS equation may be suggestive con- 



cerning the relation between the KS and KPZ equations 
in one and two spatial dimensions. Moreover, the noisy 
KS equation studied here may be relevant for the under- 
standing of dynamic roughening in, e.g., sputter eroded 
surfaces and, in principle, to any physical system mod- 
eled by the deterministic KS equation in which the rele- 
vance of time dependent noise as, e.g., fluctuations in a 
flux or thermal fluctuations, can be argued for. 

The outline of this paper is as follows. In the next sec- 
tion we introduce the noisy KS equation and discuss how 
it naturally arises in the description o f io n-sputtering 
Then we derive the RG flow in section 



III. Section 



contains the analysis of the RG flow and our results for 
one and two spatial dimensions. Finally, in section [v] we 
conclude and summarize. 



II. THE NOISY KURAMOTO-SIVASHINSKY 
EQUATION 

In this section we introduce the noisy KS equation and 
address some of its peculiarities and limiting cases. We 
discuss a physical example in which the noisy KS equa- 
tion appears naturally, namely surfaces eroded by ion- 
sputtering. 

Recently, the dynamics of surfaces undergoing ion- 
sputtering have been studied experimentally and several 
different scaling behaviors found |22]-|2~i|l . Among them 
it is notewort hy t o remark the values a = 0.2 — 0.4, and 
z = 1.6 — 1.8 J22] , which are consistent with the predic- 
tions of the KPZ equation in 2 + 1 dimensions |§||- On 
the other hand, it is well known that for amorphous tar- 
gets, ion-sputtering leads in many cases to the formation 
of a periodic ripple structure whose orientation depends 
on the angle of incidence, 9, of the ions with the normal 
to the uneroded surface |2f|. This periodic structure is 
associated with an instability in the system. An experi- 
mental study of this kind of morphologies can be found in 
p3| . In Ref. |2q] a nonlinear stochastic equation has been 
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proposed to describe the dynamics of the surface profile 
height h(x,t) of a two-dimensional surface sputtered at 
off-normal angles. The 6^0 limit (normal incidence) of 
this equation results in 



dh 



vV 2 h- K(V 2 ) 2 h+^(Vh) 2 +r)(x,t). (1) 



Here, we will consider the case where the field h(x,t) 
describes the height profile of a d-dimensional surface 
evolving in a d + f-dimcnsional medium. The surface 
tension coefficient v is negative, whereas K is a posi- 
tive surface diffusion coefficient |p7|| . The strength of 
the nonlinearity is given by A, and rj(x,t) is a Gaussian 
white noise with zero mean and short-range correlations 
described by 



(ri(x, t)r)(x',t')) = 2DS d {x - x')5{t - t'). 



(2) 



The fact that v < means that the system is linearly 
unstable, a fact which in ion-sputtered systems is related 
to faster erosion velocity at the bottom of the troughs 
than at the peaks of the crests pj|] , which in turn is re- 
lated to the formation of the periodic ripple structure re- 
ferred to above. The same kind of instability takes place 
in the deterministic KS equation, which corresponds to 
Eq. (0) without the noise term. In the following we will 
refer to Eq. (0) as the noisy KS equation. From a linear 
stability analysis one finds that the amplitude for solu- 
tions of the form 



ho(x,t) ~ exp(il/ci) sin(/cx) 



(3) 



is characterized by the rate fife = —(vk 2 + Kk ). This 
expression is plotted in Fig. Q and seen to have a positive 
value (corresponding to unstable modes in the system) 
for momenta between and fcp, where 



k 



(4) 



The modes with k > k are stable. In this sense, k 
marks the onset of the unstable modes, and the length 
scale 1/ko can be related to the wavelength of the ripple 
structure. The nonlinearity A couples the stable and un- 
stable modes, thus stabilizing the system [pT[ . If v > 0, 
it follows that all the modes in the noisy KS equation are 
stable and kg is not defined. 

If v in Eq. ([!]) is a positive coefficient, the contribution 
of the K term is expected to be negligible at large length 
scales and for long times, so that in this case the noisy 
KS equation will show the same scaling behavior as the 
KPZ equation § 



dh 



vV 2 h+ -(Vh) 2 + r)(x,t), 



(5) 



which is obtained as the K = limit of the noisy KS 
equation. However, if v < in Eq. (|l]), as is the case 
we are interested in, the surface diffusion term acts as a 



stabilizing mechanism at short length scales. In this sit- 
uation it is not a priori clear what the scaling properties 
of Eq. (§) are. 




Unstable modes 



Stable modes 



FIG. 1. Stability of Fourier modes for the noisy KS equa- 
tion. The modes from to ko are unstable, whereas those 
lying between ko and A are stable, with A the momentum cut- 
off. The shaded region corresponds to the momentum shell 
A/6 < k < A discussed in the text. 



III. RG FLOW FOR THE NOISY KS EQUATION 

The rcnormalization group is a standard tool which can 
be applied to stochastic equations in order to determine 
their scaling behavior in the large-distance long-time hy- 
drodynamic limit |^9| . Basically, the RG consists of the 
combination of a coarse-graining of the system followed 
by rescaling by a factor b = e . Successive applications 
of this transformation leads to the RG flow of the pa- 
rameters appearing in the equation in terms of the scale 
£. As usual, we will be interested in considering the vari- 
ation of the parameters when t is infinitesimal, that is, 
the flow is given by first order differential equations with 
respect to i. The critical exponents a and z characterize 
the scaling behavior of, e.g., the correlation function 

C{x -x',t- t') = ([h(x, t) - h(x', i')] 2 ) 

= \ X - X i\ 2 «g(\t-t>\/\x-x'n, (6) 

where g is a scaling function. These exponents are cal- 
culated at the fixed points of the RG transformation. 

The standard implementation of the RG procedure to 
deal with equations such as Eq. (Q) proceeds first by 
transforming it to Fourier space. In our case, it can be 
easily seen that the instability (y < 0) in Eq. (|l|) gener- 
ates a pole for zero frequency at the wavevector k = ko 
in the bare propagator 
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yk 2 + Kk 4 ' 



(7) 



As a result, when trying to solve the equation addi- 
tional divergencies arise together with those existing for 
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k — ► 0. The RG circumvents these divergencies by inte- 
grating out a small momentum shell and can consistently 
be performed. The momentum shell corresponds to the 
wavevectors 



A/6 < k < A, 



(8) 



where the momentum cutoff A = 2ir/a is related to the 
lattice spacing a, which acts as a short distance cutoff. 
Without loss of generality, wc take A = 1 in the following. 

The existence of the fco pole splits the momentum axis 
in two regions where the coarse-graining procedure can 
be applied. The cutoff A can either be put in the band 
of stable modes, cf. Fig. [I], or in the band of unstable 
modes. We will discuss the outcome of the RG analysis 
for both choices. Intuitively, however, one should coarse- 
grain over a shell of stable modes and study the effect of 
this coarse-graining on the large scales by means of the 
RG transformation. 

The calculation of the complete RG flow for the noisy 
KS equation ([!]) is tedious. We follow standard diagram- 
matic procedures |p],|3"cfl and here we only give the final 
result (see, e.g., |3f|]). Compared to the calculations for 
the KPZ equation (K = 0) [pP,|30||, the main source of 
the complications is that we have to keep terms up to or- 
der fc 4 in our expansion, in order to be able to calculate 
the renormalization of the —Kk A term. The complete 
RG flow for the noisy KS equation at one-loop order is 

dv ( A 2 L> v(2-d)+K(A-d)\ 

dK ( X 2 D 
-dl= K { Z - 4+ — K « 



clqv 3 + aiis 2 K + a 2 vK 2 + a 3 K 3 
16d(d + 2)(v + K) 5 



X(a + z - 2), 



dX 

d£ 

d l = D [ z - 2a - d+yDK U(. + Kr 



,(10) 

(11) 

(12) 



with K d = S d /{2n) d , where S d = 2n d / 2 /T{d/2) is the 
surface area of the d dimensional unit sphere. The poly- 
nomials dj = a%{d) (i = 0, 1, 2, 3) are given by 

a = 3d(d-2), oi = lid 2 - 2Ad- 20, 

a 2 = 13d 2 - AOd - 60, a 3 = 5d 2 - 22d - 16. (13) 

When studying the above RG flow, the coupling con- 
stants are conveniently taken to be 



_ IQtf_D _ K 

9v ~ Ad v 3 ' ~ v ■ 

Using these variables, the flow (p|) — (|l2|) reads 



(14) 



dg v Ad~6 + 3(d-A)f„ 2 
— = (2 - d)g v + g„ 



(15) 



dU - . b + biU + b 2 f 2 + b 3 fl + bjft 

dl ~ ZJu+ 4(d + 2)(l + /„) 8 9l " [ ' 

The polynomials bi — bi(d) (i = 0, 1, . . . , 4) are given by 

b = 3d 2 - 6d, h = lbd 2 - 2Ad - 36, 

b 2 = 25d 2 - A8d - 124, b 3 = 17d 2 - 38d - 96, 



b 4 = Ad 2 -8d- 32. 



(17) 



Note that g v is the only coupling constant one needs to 
study in the KPZ case. In our case the flow for g v is 
affected by the additional coupling /„, which probes the 
relevance at large distances of surface diffusion with re- 
spect to surface tension. The pair (f v ,g v ) is convenient 
to analyze in the cases where K is smaller than v. On the 
other hand, when v is flowing towards zero the natural 
coupling constants are 



\ 2 D 



1 



9K 



IK 



9u 



I6d(d + 2) K 3 A{d+2) /3 
K 3v ' 



(18) 



for which the RG flow becomes 



d 9K ,„ „ c + c x f K + c 2 /i + c 3 /|- 2 

-(8-%k — — ^ g K , (19) 



dl 
df K 



9f b i + b z f K + b 2 f 2 K + b 1 .f K + bof i K 

di = 2fK {iTIkJ 5 ffA '' (20) 

with the polynomials Ci = Ci(d) [i = 0, 1,2, 3) 

c = lid 2 - 7 Ad - 48, ci = 31d 2 - 136d - 180, 
c 2 = 29d 2 - 80d - 60, c 3 = 9d 2 - 18d. (21) 

The flows for Eqs. (|l||)-((l|) and Eqs. ©-© are 
singular at the lines /„ = — 1 and fx = — 1, respectively. 
This is the signature of the singularity (^J) encountered 
in the bare propagator (Q) for the finite value of the mo- 
mentum k — k . Note also that the two sets of flow 
equations describe the same system and, accordingly, we 
can analyze the effect of the RG transformation using 
any of them. 



IV. RESULTS AND DISCUSSION 

The flow equations dTJ)-© and @-@ contain 
our results for any dimension d. First, we will determine 
the fixed points (FP) for the RG flow, and the values 
of the critical exponents a and z. Then, we will con- 
sider the implications for the flow concerning the large 
distance behavior of the system when we coarse-grain 
in the different regions of momentum space depicted in 
Fig. 0. We will study the physically relevant cases of 
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one and two spatial dimensions. In Figs. 2—5, we have 
shown the RG flows obtained from the flow equations for 
d = 1,2. The behavior displayed in the figures is quali- 
tative, and the drawings are not in scale. The physical 
region of the parameter space is determined by A 2 > 0, 
D > 0, and K > 0. We consider a system which is ini- 
tially characterized by v < 0, and in order to determine 
the renormalization of the surface tension term we will 
fix the values of K and D. 



/v=-l 



a + z 




FIG. 2. Schematic RG flow for in d = 1. There 

are three fixed points: The EW fixed point at the origin, the 
KPZ attractive fixed point, and a stable spiral FP. 




FIG. 3. Schematic RG flow for (//c,5if) in d = 1. There 
are three fixed points: The linear MBE fixed point at the 
origin, the KPZ attractive fixed point, and a stable spiral FP. 



(22) 



which is valid in any dimension for any finite fixed point 
at which A ^ 0. 

In d = 1 we find three FP's, cf. Fig. |^. One is the 
unstable origin (saddle point), which is characterized by 
the Edwards- Wilkinson (EW) exponents p3| 



(23) 



corresponding to the KPZ equation with A = 0. The 
second one is a stable FP at (/*,#*) = (-0.04,0.539) [or 
(f^,g* K ) = (-25.25, -722.8)] with exponents 



2 = 1.46, a = 0.54. 



(24) 



We identify this as the KPZ fixed point and believe the 
negative value of /* (i.e., K* < 0) to result from the 
one-loop approximation. The third FP is a stable focus 
(spiral FP, characterized by eigenvalues of the linearized 
RG transformation which are imaginary numbers with 
negative real parts) at (/*,#*) = (-0.248,-1.865) [or 
g* K ) = (-4.037, 10.204)] with exponents 



z = 3.12, 



a = -1.12. 



(25) 



In d = 2 we find two fixed points, cf. Fig. [|. One is the 
unstable origin (saddle point) with the EW exponents 



2 = 2, 



a = 0. 



(26) 



The second one is also an unstable FP (saddle point) at 
(f:,9t) = (1/3,-1.757) [or (f* K ,g* K ) = (3,-2.965)] with 
exponents 



2.49, 



a = -0.49. 



(27) 



Observe that, in the (/k^k) variables a new FP ap- 
pears at the origin for any dimension, cf. Figs. || and |^, 
with the exponent values 



4-d 



(28) 



corresponding to the linear Molecular-Beam-Epitaxy 
(MBE) equation PJ-^l, which is obtained from the 
noisy KS equation (n]) by setting v = A = 0. In our case, 
this FP is always unstable, which reflects the irrelevance 
of surface diffusion at large distances in the presence of 
surface tension and a KPZ nonlinearity. 



A. Fixed Points and Critical Exponents 

As previously stated, the critical exponents are calcu- 
lated at the FP's of the RG flow (|)-((l|). Equation (0) 
reflects the symmetry of the noisy KS equation under an 
infinitesimal tilt of the interface (Galilean transforma- 
tion) Hf7],[30| and yields the exponent identity 



B. RG flow for A > k 

First, let us coarse-grain the noisy KS system in the 
band of stable modes, i.e., for A > ko. Remember that 
we use units in which A = 1, so that taking A > k is 
equivalent to taking ko < 1. In terms of the coupling 
constants and from Eq. (||) this yields \f„\ > 1, and since 



4 



v < 0, this implies — oo < f v < —1. Therefore, taking 
A > fco implies that we will take initial values of v and 
if such that 



■ oo 



< /„ < -1, g v < 0, 



(29) 



or. 



equivalently, 



-Kf K < 0, g K > 0. 



In the remaining subsections we will use the same con- 
vention, so that whenever we discuss the flow for a differ- 
ent value of A (e.g., A ^> fc ) this corresponds to taking 
initial values of v and K such that the value of fco is in 
the corresponding relation to A = 1 (e.g., fco <C 1). 



/v=-l 




FIG. 4. Schematic RG flow for (f v ,g„) in d = 2. There 
are two fixed points: The EW fixed point at the origin and a 
saddle point. The flow in the quadrant /„ , g v > is eventually 
towards the KPZ strong-coupling fixed point. 




FIG. 5. Schematic RG flow for {f K ,g K ) in d = 2. There 
are two fixed points: The linear MBE fixed point at the origin 
and a saddle point. 



1. Dimension d — 1 

We consider the flow (with initially v < 0) in the 
{fvi9v) variables, cf. Fig. ||. If we take A close to fc , 
then \v\ decreases monotonously, whereas (for a fixed D) 
the quantity |A 2 /z^ 3 increases. This flow takes us to very 
large values of /„ and g„, where the analysis would even- 
tually become inconclusive, so that the behavior is more 
adequately studied in the (Jk,9k) variables, cf. Fig. 
Before doing that, we note that the same asymptotic be- 
havior is reached if we start out with a value of A ^> fco- 
In this case, however, initially \v\ and A grow steadily, 
until A 2 reaches a value at which df v /dt — 0. Thereafter, 
the flow merges with the one already considered for an 
initial value A > fco- 

Looking at the (fa, 9k) variables, we find that starting 
with A > fco, the surface tension increases monotonously 
(i.e., it becomes less negative), whereas A 2 decreases till it 
reaches a line on which dgx/di = 0. From that moment 
on, both A 2 and v increase, eventually crossing the v = 
axis (Jk — 0) and yielding a flow towards the KPZ fixed 
point. If we take A 3> fco we reach similar conclusions. 
First, \v\ and A 2 increase until the flow turns back on a 
line along which d/x/di = 0. Thereafter, the flow merges 
with that studied for A > fco, and crosses to positive v 
values. Now, both fx and gx grow steadily, so that 
relevant information is gained about later stages of the 
flow if we come back to the (/ y ,<7^) variables. In the 
quadrant f u ,9v > for d = 1 (cf. Fig. ^), we observe 
that the flow is towards the fixed point characterized by 
the KPZ exponents, cf. Eq. (§|). 

To summarize, the flow which initially started out with 
a negative value of v has renormalized towards the KPZ 
fixed point, which is characterized by a positive value 
of v, when we take the momentum cutoff to lie in the 
band of stable modes and coarse-grain the system in that 
region. This behavior is in qualitative agreement with the 
results reported in Refs. JT2| , pT[ -|l^] , which investigated 
the deterministic KS equation by various coarse-graining 
procedures different from the method used in the present 
paper. Also, it has been argued in Ref. |[?J that for a 
restricted range of initial parameter values the noisy KS 
equation scales to a KPZ equation, which is in accordance 
with our results obtained by investigating the complete 
set of flow equations. 

It is instructive to study the evolution under the RG 
flow of the finite pole 



(30) 



associated with the instability in the system, cf. the dis- 
cussion following Eq. (Q). A naive scaling argument 
(which would be correct for the linear equation) leads 
to a transformation under rescaling with a factor b = e £ 
as 




K 



(31) 
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implying that under successive rescalings one would get 
kg — > bko — > oo, and we would be left with unstable 
modes only. However, we have seen that under the RG 
flow for the nonlinear equation, at some point the surface 
tension v becomes zero. As a result, fco — * 0, so that 
under the RG flow the band of stable modes shrinks to 
zero, and the noisy KS system evolves as the stable KPZ 
equation. 



2. Dimension d - 



2. Dimension d = 2 

For a two-dimensional interface, we can see in Fig. 5| 
that the flow starting in the region ( |32| ) points towards 
the origin, which as noted above (see Eq. (^)) is char- 
acterized by the EW exponents in 2+1 dimensions, i.e., 
the interface is flat. Interestingly, this is a similar scaling 
behavior to that found in Ref. [^0| for the scaling solution 
of the deterministic KS equation in 2+1 dimensions, in 
that the same value z = 2 for the dynamic exponent is 
also obtained. 



The analysis of the RG flow for A > ko in d = 2 leads to 
analogous conclusions to the one-dimensional case. We 
study the (f v ,9v) variables for — oo < /„ < —1. The flow 
in this case is completely analogous to the d = 1 case 
(compare the corresponding regions in Figs. || and so 
that {fui9i>) renormalize to very large negative values. 
Again, we switch to the {Jk ,9k) flow, cf. Fig. ||, where 
we observe that v is crossing the zero axis, thus renormal- 
izing to a positive value. Once we are in the fx > 9k > 
region, the flow is towards large positive values for fx, 
gx, and we return to the f v ,9v > variables. Here, we 
observe that f v flows towards the f v = axis (which 
cannot be crossed), while g v steadily increases. This we 
interpret as the irrelevance of surface diffusion at large 
distances, where the flow is eventually towards the KPZ 
strong coupling fixed point, unaccessible to the dynamic 
RG perturbative approach 0. 



C. RG flow for A < k 

We will now analyze the behavior of the noisy KS sys- 
tem when we integrate out a momentum shell in the 
band of unstable modes. The regions in the (fv,gv) and 
{.Ik , gx) planes which correspond to taking a value for 
the momentum cutoff lying in the band of unstable modes 
are 



- 1 < /„ < 0, g v < 0, 
-oo < fx < -1, gx > 0, 



(32) 



cf. the discussion in the beginning of subsection IV B| . 

Taking A < fco means that all the modes in the system 
are linearly unstable. Consequently, it would be desirable 
to check the results reported in this subsection through 
some technique of a different nature, as, e.g., a numerical 
integration of the noisy KS equation. 



1. Dimension d — 1 

In this case, as can be seen in Fig. |[ the flow is at- 
tracted by the stable focus. The negative roughness ex- 
ponent a that characterizes this FP (cf. Eq. (p5[)) could 
indicate a flat interface. 



D. RG flow for the region v > 0, K < 

For the sake of completeness, we devote the last two 
subsections to the analysis of the RG flow in the other 
parameter regions of Figs. ||— [|. First, we will consider 
the case in which the initial values of the parameters are 
such that v > and K < 0. In terms of the initial values 
of the coupling constants, this corresponds to the region 



f v < 0, g u > 0, 
fx < 0, g K < 0. 



(33) 



In this case, the band of stable modes is [0, fco], whereas 
the unstable modes are those with k > fco, that is, the 
stability of the two momentum regions has been inter- 
changed with respect to the discussion in subsections 
fVB| and |IV C| . Therefore, if we consider A < fco, this 
means that —1 < f u < 0, and, moreover, that A lies in 
the band of stable modes. 

In one dimension, when A ^ fco the flow is towards the 
singular line f v — — 1 , whereas for A <C fco the flow can in 
principle reach the KPZ fixed point, cf. Fig. [| The latter 
is the behavior one expects for a KPZ equation in which 
one has included a negative surface diffusion coefficient, 
since such a term would be irrelevant with respect to the 
stable surface tension term. Analogous conclusions can 
be drawn in the d = 2 case (cf. Fig. I|), for which taking 
A < fco again leads to a flow towards the line /„ = — 1, 
whereas for A <C fco the flow is towards a large positive 
value of <?„, which we interpret as pointing towards the 
strong coupling KPZ behavior. 

On the other hand, if we coarse grain a system for 
which A > fco, we place ourselves in the band of unstable 
modes. Both in one and two dimensions, it can be seen 
that fco flows to a finite value, 



1 (=A), 



(34) 



as follows, e.g., from Figs. ^| and |] in the region where 
— oo < f v < — 1, g v > 0. Consequently, the singularity 
remains at a fixed value of the momentum, and it is 
reached at some stage along the coarse-graining proce- 
dure. In the flow diagrams, it can be observed that the 
flow terminates on the singular line /„ = fx = — 1, both 
for d = 1, 2, where the analysis is inconclusive. 



G 



E. RG flow for the region /„ > 0, g v < 

Finally, we discuss the flow in the region given by 
f u > 0, g u < 0, which corresponds to the following pos- 
sibilities for the parameters appearing in the noisy KS 
equation: i) If both v, K > 0, then A 2 < 0, or D < 0, 
which is unphysical. ii) If both v, K < 0, this means we 
are dealing with a highly linearly unstable equation. 

In d — 1, fixing K the flow is initially towards small 
values of /„ (i.e., v — > — oo), while A 2 increases steadily. 
The later flow is towards large values of /„, g v , so we 
study it in the {fK,gx) variables. In terms of these (cf. 
Fig. ||), we observe that for fx ~ 0, gx < the flow is 
ultimately towards the singular line Jk = — 1, where the 
analysis is inconclusive. 

For d = 2, there exist two possibilities separated by the 
saddle FP, cf. Fig. ||. If the starting {f v ,9v) point lies to 
the right of the stable separatrix, the flow is eventually 
towards the fx = — 1 singular line, following a similar 
reasoning to that in the d — 1 case. However, if the 
initial value of {f u ,g w ) is above the separatrix (as, e.g., 
starting out with a small f v , corresponding to a huge 
value of \v\ or to a very small value of K), the flow is 
towards the origin, characterized by Edwards- Wilkinson 
exponents in 2+1 dimensions. 



V. CONCLUSIONS 

In the present paper we have introduced a noisy version 
of the Kuramoto-Sivashinsky equation, which we have in- 
vestigated by means of a dynamic renormalization group 
procedure. We have analyzed simultaneously two sets of 
variables which are combinations of the parameters ap- 
pearing in the equation. For a system in which the lattice 
spacing is smaller than the typical wavelength of the lin- 
ear instability occurring in the system, the large-distance 
and long-time behavior of this equation is the same as for 
the KPZ equation in one and two spatial dimensions. For 
the d = 2 case the agreement is only qualitative, since due 
to the strong coupling behavior there does not exist a (fi- 
nite) fixed point accessible to the dynamic RG perturba- 
tive approach. On the other hand, when coarse-graining 
on larger scales the asymptotic flow depends on the initial 
values of the parameters, and no general conclusions can 
be obtained. It would be interesting to check the results 
obtained here, especially those for A < ko, where all the 
modes in the system are unstable, with the outcome of a 
qualitatively different approach such as, e.g., a numerical 
integration of the noisy Kuramoto-Sivashinsky equation. 

Finally, given the fact that the noisy KS equation con- 
tains the main physical mechanisms relevant to the dy- 
namics of sputter-eroded surfaces, we believe that the 
results obtained here may be relevant for the description 
of the roughening of such surfaces. Note, however, that 
the noisy KS equation is the isotropic limit (9 = 0) of 
the equation derived in Ref. |2q]. The consequence of 



this limit for the large distance behavior of the system is 
a question that should be elucidated by the investigation 
of the full equation proposed in Ref. p6j, of which the 
present work can be thought of as a first step. 
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